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Abstract 

The inclusive charged particle transverse momentum distribution is measured in proton-proton collisions at y[s = 
900 GeV at the LHC using the ALICE detector. The measurement is performed in the central pseudorapidity region 
(I77I < 0.8) over the transverse momentum range 0.15 < pj < 10 GeV/c. The correlation between transverse mo- 
mentum and particle multiplicity is also studied. Results are presented for inelastic (INEL) and non-single-diffractive 
(NSD) events. The average transverse momentum for < 0.8 is {p T } INE l = 0.483+0.001 (stat.) +0.007 (syst.) GeV/c 
and (p T ) nsd = 0.489 + 0.001 (stat.) +0.007 (syst.) GeV/c, respectively. The data exhibit a slightly larger (p T ) than 
measurements in wider pseudorapidity intervals. The results are compared to simulations with the Monte Carlo event 
generators PYTHIA and PHOJET. 
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1. Introduction 

The precise measurement of the transverse momen- 
tum spectrum of charged particles produced in proton 
collisions in the energy range of the Large Hadron Col- 
lider (LHC) HI offers unique information about soft 
and hard interactions. Perturbative Quantum Chromo- 
dynamics (pQCD) is a framework for the quantitative 
description of parton-parton interactions at large mo- 
mentum transfers, i.e. hard scattering processes. How- 
ever, a significant fraction of the particles produced in 
pp collisions do not originate from hard interactions, 
even at LHC energies. In contrast to hard processes, 
the description of particle production in soft interactions 
is not well-established within QCD. Current models of 



hadron-hadron collisions at high energies, such as the 
event generators PYTHIA 0] and PHOJET [3], com- 
bine perturbative QCD for the description of hard par- 
ton interactions with phenomenological approaches to 
model the soft component of the produced particle spec- 
trum. Data on charged particle production in hadron- 
hadron collisions will have to be used to tune these mod- 
els before they can provide a detailed description of the 
existing measurements and predictions for particle pro- 
duction characteristics in pp collisions at the highest 
LHC energies. These data include the measurement of 
multiplicity, pseudorapidity (rj) and transverse momen- 
tum (pr) distributions of charged particles and correla- 
tions, such as the dependence of the average transverse 
momentum, (pr), on the charged particle multiplicity. 
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The charged particle pseudorapidity densities and 
multiplicity distributions in pp collisions at -\fs = 0.9, 
2.36 and 7 TeV were presented in recent publications 
by the ALICE collaboration J3, S Hi- In this letter, 
we present a measurement in pp collisions at -\fs = 
900 GeV of the transverse momentum spectrum of 
primary charged particles and the correlation between 
(pr) and the charged particle multiplicity. Primary 
particles include particles produced in the collision or 
their decay products, except those from weak decays of 
strange hadrons. The measurement is performed in the 
central rapidity region (|^| < 0.8) and covers a pj range 
0.15 < pj < 10 GeV/c, where both hard and soft pro- 
cesses are expected to contribute to particle production. 
The data from the ALICE experiment presented in this 
letter serve as a baseline for future studies of pp colli- 
sions at higher LHC energies and particle production in 
heavy-ion collisions [0. 

2. Experiment and data collection 

The data were collected with the ALICE detec- 
tor Jit] during the startup phase of the LHC in December 
2009. The ALICE detector, designed to cope with high 
track densities in heavy-ion collisions, provides excel- 
lent track reconstruction and particle identification ca- 
pabilities. This also makes the detector well-suited to 
detailed studies of global characteristics of pp interac- 
tions |0. 

In this analysis of the first pp collisions at -\fs = 
900 GeV, charged particle tracking and momentum re- 
construction are based on data recorded with the Time 
Projection Chamber (TPC) and the Inner Tracking Sys- 
tem (ITS), both located in the central barrel of ALICE. 
The detectors in the central barrel are operated inside a 
large solenoidal magnet providing a uniform 0.5 T field. 

The ALICE TPC H is a large cylindrical drift de- 
tector with a central high voltage membrane maintained 
at -100 kV and two readout planes at the end-caps. 
The active volume is limited to 85 < r < 247 cm and 
-250 < z < 250 cm in the radial and longitudinal 
directions respectively. The material budget between 
the interaction point and the active volume of the TPC 
corresponds to 11% of a radiation length, averaged in 
\r]\ < 0.8. The central membrane at z = divides the 
nearly 90 m 3 active volume into two halves. The ho- 
mogeneous drift field of 400 V/cm in the Ne-C02-N 2 
(85.7%-9.5%-4.8%) gas mixture leads to a maximum 
drift time of 94 ps. Ionization electrons produced by 
charged particles traversing the TPC drift towards the 
readout end-caps composed of 72 multi-wire propor- 
tional chambers with cathode pad readout. The typi- 



cal gas gain is 10 . Signals induced on the segmented 
cathode planes, comprising a total of 558k readout pads, 
are transformed into differential semi-gaussian signals 
by a charge-sensitive shaping amplifier (PASA). This is 
followed by the ALICE TPC ReadOut (ALTRO) chip, 
which employs a 10 bit ADC at 10 MHz sampling rate 
and four digital filtering circuits. These filters also per- 
form tail cancellation and baseline restoration. They are 
optimized for precise position and dE/dx measurements 
in the high track density environment of heavy-ion colli- 
sions. To ensure optimal drift and charge transport prop- 
erties, the TPC was operated with an overall tempera- 
ture uniformity of AT m 60 mK (r.m.s.). The oxygen 
contamination was less than 5 ppm. 

The ITS is composed of high resolution silicon 
tracking detectors, arranged in six cylindrical layers at 
radial distances to the beam line from 3.9 to 43 cm. 
Three different technologies are employed. 

For the two innermost layers Silicon Pixel Detec- 
tors (SPD) are used, covering the pseudorapidity ranges 
\rj\ < 2 and \tj\ < 1.4, respectively. A total of 9.8 mil- 
lion 50 x 425 pm 2 pixels enable the reconstruction of 
the primary event vertex and the track impact parame- 
ters with high precision. The SPD was also included in 
the trigger scheme for data collection. 

The SPD is followed by two Silicon Drift Detec- 
tor (SDD) layers with a total of 133k readout channels, 
sampling the drift time information at a frequency of 
20 MHz. The SDD are operated with a drift field of 
500 V/cm, resulting in a drift speed of about 6.5 pm/ns 
and in a maximum drift time of about 5.3 ps. 

The two outermost Silicon Strip Detector (SSD) lay- 
ers consist of double-sided silicon micro-strip sensors 
with 95 pm pitch, comprising a total of 2.6 million read- 
out channels. Strips of the two sensor sides form a 
stereo angle of 35 mrad, providing two-dimensional hit 
reconstruction. 

The design spatial resolutions of the ITS sub- 
detectors {cr r< p x 0%) are: 12 x 100 pm 1 for SPD, 
35 x 25 pm 2 for SDD, and 20 x 830 pm 2 for SSD. The 
SPD and SSD detectors were aligned using survey mea- 
surements, cosmic muon data [ 10] and collision data to 
an estimated accuracy of 10 pm for the SPD and 15 pm 
for the SSD. No alignment corrections are applied to the 
positions of the SDD modules, for which calibration and 
alignment are in progress. The estimated misalignment 
of the SDD modules is about 100 pm. The TPC and 
ITS are aligned relative to each other to the level of a 
few hundred micrometers using cosmic-ray and pp data 
by comparing pairs of track segments independently re- 
constructed in the two detectors. 

The two forward scintillator hodoscopes (VZERO) 
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are included in the trigger. Each detector is segmented 
into 32 scintillator counters which are arranged in four 
rings around the beam pipe. They are located at dis- 
tances z — 3.3 m and z = -0.9 m from the nominal 
interaction point and cover the pseudorapidity ranges: 
2.8 < rj < 5.1 and -3.7 < 77 < -1.7 respectively. The 
time resolution of about 1 ns of the VZERO hodoscope 
also allows for a discrimination against beam-gas inter- 
actions. 

During the startup phase of the LHC in 2009, 
four proton bunches per beam were circulating in the 
LHC with two pairs of bunches crossing at the AL- 
ICE intersection region and protons colliding at -\fs = 
900 GeV. The detector readout was triggered using the 
LHC bunch-crossing signals in coincidence with signals 
from the two upstream beam pick-up counters and a 
minimum-bias interaction trigger requiring a signal in at 
least one of the SPD pixels or one of the VZERO coun- 
ters J3, S . Events with only one bunch or no bunches 
passing through ALICE were also recorded to study 
beam related and random background. 

3. Data analysis 

The total inelastic pp cross section is commonly 
subdivided into contributions from diffractive and non- 
diffractive processes. To facilitate comparison with 
existing measurements, we perform our analysis for 
two classes of events: inelastic (INEL) and non-single- 
diffractive (NSD) pp collisions. 

In this analysis, 3.44 x 10 s triggered pp events at 

= 900 GeV are analyzed. To remove beam related 
background events, an offline event selection based on 
the VZERO timing signal and the correlation between 
the number of hits and tracklets in the SPD is applied as 
in 0], reducing the sample to 2.67 X 10 5 events. This 
event selection is refered to as MBor JH. 

For the INEL analysis we use the event sample se- 
lected with the MBor condition. A subset of these 
events (2.15 x 10 5 ) is used for the NSD analysis, se- 
lected offline by requiring a coincidence between the 
two VZERO detectors (the MBand selection). This 
condition suppresses a significant fraction of the single- 
diffractive events and hence reduces the systematic er- 
rors related to model dependent corrections ||5J ■ 

The fractions of the different process types con- 
tributing to the selected event samples are estimated by 
a Monte Carlo simulation, implementing a description 
of the ALICE detector response ll 111 to pp collisions 
at V* = 900 GeV from the PYTHIA event generator 
version 6.4.21 tune D6T (109) ll 211 . The process frac- 
tions of single-diffractive (SD) and double-diffractive 



(DD) events are scaled in Monte Carlo to match the 
cross sections in pp at yfs — 900 GeV measured by 
UA5 ifijll . The selection efficiency for INEL events us- 
ing MBor and NSD events using MBand is approxi- 
mately 96% and 93%, respectively @|- 

Charged particle tracks are reconstructed using in- 
formation from the TPC and ITS detector systems. Sig- 
nals on adjacent pads in the TPC are connected to parti- 
cle tracks by employing a Kalman filter algorithm. The 
TPC tracks are extrapolated to the ITS and matching hits 
in the ITS detector layers are assigned to the track. In 
order to maximize the hit matching efficiency and avoid 
possible biases of the track parameters due to the non- 
uniform degree of alignment of the ITS sub-detectors, 
the space point uncertainties of the ITS hits are set to 
100 fxm for SPD and 1 mm for both SDD and SSD. 

The event vertex is reconstructed using the com- 
bined track information from TPC and ITS. The tracks 
are extrapolated to the intersection region and the po- 
sition of the event vertex is fitted, using the measured 
average intersection profile as a constraint. The profile 
of the intersection region is determined on a run-by-run 
basis in a first pass through the data using the mean and 
the spread of the distribution of the reconstructed ver- 
tices. The event vertex distribution is found to be Gaus- 
sian with standard deviations of approximately 210 fim, 
250 pm, and 4.1 cm, along x, y (transverse to the beam- 
axis) and z respectively. For events where only one track 
is found, the vertex is determined from the point of clos- 
est approach of the track to the beam axis. If no track 
is found in the TPC, the event vertex reconstruction is 
based on tracklets built by associating pairs of hits of the 
two innermost ITS layers (SPD). An event with a recon- 
structed vertex position z v is accepted if |z v — Zol < 10 cm, 
corresponding to about 2.5 standard deviations of the re- 
constructed event vertex distribution centered at zo JH. 

The vertex position resolution depends on the 
event multiplicity. It can be parametrized as 
540 pm/(N sm ) 0A5 in x and y, and 550 pm/iNsm) 06 in 
z, where Nsm corresponds to the number of SPD track- 
lets. This resolution is consistent with Monte Carlo sim- 
ulations. The probability of multiple interactions in the 
same bunch crossing (pile-up) in the present data set is 
10 -4 and therefore neglected. 

The fraction of selected events in the 
MBor (MBand ) sample where an event vertex is suc- 
cessfully reconstructed is 80% (92%), resulting in a 
sample of 2.13 x 10 5 INEL (1.98 x 10 5 NSD) events 
used in the present analysis. Events, where no vertex 
is found, are included when normalizing the results. 
In order to understand and to subtract possible beam- 
induced background, the detector was also triggered 
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on bunches coming from either side of the interaction 
region, but not colliding with another bunch. From 
the study of these events we estimate that 21% of the 
triggered MBor and MBand events without a recon- 
structed event vertex or with zero selected tracks are 
background, and the number of events used for normal- 
ization of the final results is corrected accordingly. The 
estimated contribution from beam-induced background 
events to the event sample, where a vertex was found, is 
negligible. From the analysis of empty bunch crossing 
events the random contribution from cosmics and noise 
triggers is also found to be negligible. 

To study the transverse momentum spectrum, 
charged particle tracks are selected in the pseudorapid- 
ity range \rj\ < 0.8. In this range, tracks in the TPC 
can be reconstructed with maximal length, and there 
are minimal efficiency losses due to detector bound- 
aries. Additional quality requirements are applied to 
ensure high tracking resolution and low secondary and 
fake track contamination. A track is accepted if it has 
at least 70 out of the maximum of 159 space points in 
the TPC, and the x 1 P er space point used for the mo- 
mentum fit is less than 4. Additionally, at least two 
hits in the ITS must be associated with the track, and 
at least one has to be in either of the two innermost lay- 
ers, i.e., in the SPD. The average number of associated 
hits per track in the six ITS layers is 4.7, mainly deter- 
mined by the fraction of inactive channels, and is well 
reproduced in Monte Carlo simulations. Tracks with 
Pt < 0.15 GeV/c are excluded because their recon- 
struction efficiency drops below 50%. Tracks are also 
rejected as not associated to the primary vertex if their 
distance of closest approach to the reconstructed event 
vertex in the plane perpendicular to the beam axis, do, 
satisfies do > 0.35 mm + 0.42 mm x p T 09 , with pj in 
GeV/c. This cut corresponds to about seven standard 
deviations of the pj dependent transverse impact pa- 
rameter resolution for primary tracks passing the above 
selection. It is tuned to select primary charged particles 
with high efficiency and to minimize the contributions 
from weak decays, conversions and secondary hadronic 
interactions in the detector material. The accepted num- 
ber of charged particles per event which fulfill these 
conditions is called « acc . 

With this selection, the reconstruction efficiency for 
primary charged particles and the remaining contamina- 
tion from secondaries as a function of pj are estimated 
by Monte Carlo simulation using PYTHIA, combined 
with detector simulation and event reconstruction. The 
procedure estimates losses due to tracking inefficiency, 
charged particles escaping detection due to weak de- 
cay, absorption and secondary interaction in the detec- 



tor. The inefficiencies of the event selection and of the 
event vertex reconstruction are accounted for. The lat- 
ter two affect mostly low-multiplicity events, which im- 
poses a bias on the uncorrected pj spectrum due to the 
correlation between multiplicity and average momen- 
tum. 
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Figure 1 : Charged-particle track reconstruction efficiency for primary 
particles (top) and contamination from secondary particles (bottom), 
for positively and negatively charged particles in \t]\ < 0.8 as a func- 
tion of pt- The tracking efficiency is normalized to the number of 
generated primary particles using PYTHIA. The contamination from 
secondary tracks was scaled in Monte Carlo to match the measured do 
distributions (see text). 

The primary charged particle track reconstruction 
efficiency in the region \t]\ < 0.8 reaches 75% at pj ~ 
1 GeV/c, as shown in Fig.Q] The slight decrease of effi- 
ciency observed for pj > 1 .5 GeV/c is a consequence of 
the projective segmentation of the readout plane in az- 
imuth, causing stiff tracks to remain undetected if they 
fall between two adjacent TPC readout sectors. For 
pj < 0.6 GeV/c, the reconstruction efficiency decreases 
and reaches 50% at 0.15 GeV/c. The losses at low 
Pt are mainly due to energy loss in the detector ma- 
terial and to the track bending in the magnetic field. No 
significant dependence of the track reconstruction effi- 
ciency on the track density is observed in simulations 
for charged particle multiplicities relevant for this anal- 
ysis. The contamination from secondary particles such 



9 



as charged particles from weak decays, electrons from 
photon conversions, and products from secondary inter- 
actions in the detector material is also shown in Fig.Q] 
It has a maximum of 9% at the lowest pj and drops 
below 3% for pj > 1 GeV/c. A comparison of the 
c/o distributions of data and Monte Carlo tracks indi- 
cates that the Monte Carlo simulation using PYTHIA 
underestimates the particle yield from secondaries by 
0-50%, depending on pj. This is consistent with the 
fact that PYTHIA underestimates the strangeness yield 
by a similar amount, when compared to previous results 
in pp and pp collisions HUH]. For the final correc- 
tions to the data we scale accordingly the contamination 
level obtained with PYTHIA, resulting in an additional 
0-1.5% decrease of the primary particle yields. The un- 
certainty in the strangeness yield is taken into account 
in the evaluation of the overall systematic uncertainties, 
as discussed below. 

The reconstruction efficiency and contamination are 
converted to pj dependent correction factors used to 
correct the raw pj spectrum. We note that efficiency 
and secondary contamination are slightly different for 
positively and negatively charged particles, mainly due 
to the larger absorption of negatively charged particles 
and isospin effects in secondary interactions. 

The charged particle transverse momenta are mea- 
sured in the TPC, taking into account energy loss based 
on the PID hypothesis from TPC dE/dx and the mate- 
rial budget in front of the TPC. The material budget is 
studied via the measurement of electron-positron pairs 
in the TPC from photon conversions. The radial distri- 
bution of the reconstructed photon conversion points is 
compared to Monte Carlo simulations. The sum of all 
positive and negative deviations is +4.7% and -7.2%, 
respectively. The remaining material budget uncertainty 
enters into the final systematic uncertainties. In this 
analysis, we use the measurement of the momentum at 
the event vertex. 

At the present level of calibration, the transverse 
momentum resolution achieved in the TPC is given by 
(o-(p T )/p T ) 2 = (0.01) 2 +(0.007-p r ) 2 , with p T inGeV/c. 
The transverse momentum resolution for pj > 1 GeV/c 
is measured in cosmic muon events by comparing the 
muon momenta reconstructed in the upper and lower 
halves of the TPC. For p T < 1 GeV/c, the Monte Carlo 
estimate of ct(pt)/pt ~ 1% is cross-checked using 
the measured invariant mass distribution. A Monte 
Carlo based correction is applied to the pj spectra to ac- 
count for the finite momentum resolution. The correc- 
tion increases with pj and reaches 1.2% at 10 GeV/c. 

The calibration of the absolute momentum scale is 
verified employing the invariant mass spectra of A, A, 



and <f>. The reconstructed peak positions agree with 
their PDG values within 0.3 MeV/c 2 . As a cross-check, 
the qlpr distributions of particles with charge q in data 
and Monte Carlo simulation are compared and the sym- 
metry of the minimum around q/pr = is studied. 
Based on these studies, we estimate an upper limit on 
the systematic uncertainty of the momentum scale of 
\MPt)/Pt\ < 0.003. Within the p T reach of this study, 
the effect of the momentum scale uncertainty on the fi- 
nal spectra is found to be negligible. 

For the normalization of the transverse momentum 
spectra to the number of events, multiplicity depen- 
dent correction factors are derived from the event selec- 
tion and vertex reconstruction efficiencies for INEL and 
NSD events, evaluated with the PYTHIA Monte Carlo 
event generator. 

The fully corrected pjspectra are fitted by the mod- 
ified Hagedorn function 11611 

1 d2j Vch PT ( ^ PT Y b ... 

inpj dr\ dpj nij \ pro / 

For the transverse mass mj = ■Jrr% + p^. , the pion 
mass is assumed for all tracks. At small pj, the term 

v + pro) b enaves like an exponential in pj with in- 
verse slope parameter pT,o/b. This provides a good de- 
scription of the soft part of the spectrum, allowing for 
an extrapolation of the measured data to pj = 0. To as- 
sess the tail of the spectrum at pj > 3 GeV/c, a power 
law fit is performed 



2npT drj dpj 

yielding a very good description of the hard part of the 
spectrum characterized by the power n. 

The calculation of (pj) in all INEL and NSD events 
is performed using the weighted average over the mea- 
sured points in the range 0.15 < pj < 10 GeV/c com- 
bined with the result of the Hagedorn fit to extrapolate 
to pj = 0. 

In order to analyze the behaviour of (pj) as function 
of multiplicity, the INEL data sample is subdivided into 
bins of « acc . The results for (pr) are presented calculat- 
ing the weighted average over two different pj ranges, 
0.15 < p T < 4 GeV/c and 0.5 < p T < 4 GeV/c. In addi- 
tion, results are presented employing the extrapolation 
to pj — as described above. 

To extract the correlation between (pj) and the 
number of primary charged particles (n c h) in I77I < 0.8, 
the following weighting procedure is applied to account 
for the experimental resolution of the measured event 
multiplicities: 
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(Pt) («ch) = ^ (Pt) («acc)fl("ch, n. dcc ). (3) 

This method employs the normalized response ma- 
trix R(n c h, « acc ) from Monte Carlo simulations which 
contains the probability that an event with multiplicity 
« c h is reconstructed with multiplicity « acc . The results 
from this approach are consistent with an alternative 
Monte Carlo based procedure, where an average mul- 
tiplicity (« c h) is assigned to every measured multiplicty 



4. Systematic uncertainties 

In order to estimate the systematic uncertainties of 
the final pj spectra, the results of the data analysis and 
of the evaluation of the corrections from Monte Carlo 
simulations are checked for stability under varying cuts 
and Monte Carlo assumptions, within reasonable lim- 
its. In particular we studied a variation of the ratios of 
the most abundant primary charged particles (p, n, K) 
by ±30% with respect to their PYTHIA values, the rel- 
ative fractions of diffractive processes corresponding to 
their experimental errors j5[[T3t], the TPC readout cham- 
ber alignment (±100 pm), and track and event quality 
cuts in the analysis procedure. Particular attention was 
paid to the rejection efficiency of secondary particles us- 
ing the do cut. The stability of the results under vari- 
ation of the do cut value (±3 standard deviations with 
respect to the nominal value), the secondary yield from 
strange hadron decays (±30%) and the material budget 
(±10%) was studied and the systematic uncertainty is 
estimated accordingly. Systematic uncertainties of the 
ITS and TPC detector efficiencies are estimated by a 
comparison of the experimental ITS-TPC track match- 
ing efficiency with the Monte Carlo one. The systematic 
uncertainty of the VZERO triggering efficiency is stud- 
ied by varying the calibration and threshold settings in 
the data and in the Monte Carlo simulation. The event 
generator dependence is determined from a comparison 
of the PYTHIA results with those obtained using PHO- 
JET. The total systematic uncertainty on the pj spec- 
tra derived from this study is 3.0-7.1% for INEL events 
and 3.5-7.2% for NSD events, in the pj range from 
0.2 - 10 GeV/c (see Tabled). 

Also listed in Table Q] are the systematic errors in 
(pr) arising from these contributions. We note that 
only pr dependent errors on the pj spectra contribute 
to the systematic error in (pr)- Additional systematic 
uncertainties in {pj} arise from the specific choice of 



the fit function used for the pj — extrapolation, and 
the weighting procedure which is employed to derive 
(pr) as function of n c h. To estimate the uncertainty in 
the extrapolation to pj — the results are compared 
to those obtained from a fit of the Tsallis function §Y~k , 
or by fitting the spectral shape predicted by PYTHIA 
and PHOJET to our low pj data points. Based on this 
comparison a systematic error of 1 % in (py) is assigned 
to the pj = extrapolation. The weighting procedure 
(Eq. |3]l was studied using PYTHIA and PHOJET simu- 
lations. For both models, the true (pr) dependence on 
« c h from Monte Carlo can be recovered within 3% from 
the reconstructed dependence of {pj) on « acc using (f3). 
No significant multiplicity dependence of the system- 
atic errors on {pj) is observed. The total systematic 
uncertainties on (pr) are listed in Table [TJ 

5. Results and discussion 

The normalized differential yield in INEL pp colli- 
sions at tJs = 900 GeV and the fit with the parametriza- 
tion given in Eq. Q] are shown in Fig. [2] The mod- 
ified Hagedorn fit provides a good description of the 
data for pj < 4 GeV/c. The fit parameters for INEL 
events are p T ,o = 1.05 ± 0.01 (stat.)±0.05 (syst.) GeV/c 
and b = 7.92 ± 0.03 (stat.)±0.02 (syst.). The 
average transverse momentum including the extrap- 
olation to pr = is (pr) inel = 0.483 ± 
0.001 (stat.) ±0.007(syst.) GeV/c. For NSD events we 
obtain p T ,o = 1.05 ± 0.01 (stat.)±0.05 (syst.) GeV/c, 
b = 7.84 ± 0.03 (stat.)±0.02 (syst.) and (p T ) NSD = 
0.489 ± 0.001 (stat.) ±0.007 (syst.) GeV/c. Restriction 
of the modified Hagedorn fit to pj < 4 GeV/c has a neg- 
ligible effect on these results. Fig.[2]also shows the re- 
sult of a power law fit (Eq.[2) to the INEL data for pj > 
3 GeV/c. The power law fit provides a significantly bet- 
ter description of the high pj- tail of the spectrum than 
the modified Hagedorn parametrization. The result of 
the power law fit is n = 6.63 ± 0.12 (stat.) ±0.01 (syst.) 
for both INEL and NSD events. The power law shape 
of the high pj part of the spectrum is suggestive of 
pQCD. Estimates of differential cross sections can be 
obtained using the cross sections derived from the mea- 
surement by UA5 lfl3ll in pp at ^fs — 900 GeV, cr inel = 
50.3 ±0.4 (stat.) ± 1 (syst.) mb and o"nsd — 

42.6±1.4mb 

(see also IU8I0 . 

The transverse momentum distribution for NSD 
events is shown in Fig.|3](left panel) together with data 
recently published by ATLAS lH] and CMS Q, mea- 
sured in larger pseudorapidity intervals. Below pj = 
1 GeV/c the data agree. At higher pt the data are 
slightly above the other two LHC measurements. The 
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Figure 2: Normalized differential primary charged particle yield in INEL pp collisions at yfs = 900 GeV, averaged in \rf\ < 0.8. The fit ranges are 
0.15 < pr < 10 GeV/c for the modified Hagedorn function (Eq.[T} and 3 < pj < 10 GeV/c for the power law (Eq.|5J. In the lower panels, the 
ratios fit over data are shown. The open symbols indicate here data points which are not included in the fit. Errors bars are statistical only. Indicated 
as shaded areas are the relative systematic data errors. 
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Figure 3: Left panel: Normalized differential primary charged particle yield in NSD pp collisions at V* = 900 GeV, averaged in \t]\ < 0.8. The 
ALICE data are compared to results from ATLAS and CMS in pp at the same energy I19ll2rj|l . Right panel: Normalized invariant primary charged 
particle yield in NSD pp collisions at V* = 900 GeV, averaged in \rj\ < 0.8. The ALICE data are compared to results from UA1 in pp at the same 
energy 12111 . For the computation of the invariant yield, it has been assumed that all particles are pions. The shaded areas indicate the statistical and 
systematic errors added in quadrature. 
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observation of a harder spectrum is related to the differ- 
ent pseudorapidity windows (see below). 

In the right panel of Fig. [3j the normalized invari- 
ant yield in NSD events is compared to measurements 
of the UA1 collaboration in pp at the same energy 112 111 , 
scaled by their measured NSD cross section of 43.5 mb. 
As in the previous comparison to ATLAS and CMS, the 
higher yield at large pj may be related to the different 
pseudorapidity acceptances. The excess of the UA1 data 
of about 20% at low pj is possibly due to the UA1 trig- 
ger condition, which suppresses events with very low 
multiplicity, as pointed out in Il9ll . 

The results for (pr) in INEL and NSD events are 
compared to other experiments 

Em EM HQ 

in Fig. [4] Our results are somewhat higher than pre- 
vious measurements in pp and pp at the same energy, 
but in larger pseudorapidity windows. This is consistent 
with the comparison of the spectra in Fig. [3] A similar 
trend exhibiting a larger (pr) in a smaller pseudorapid- 
ity interval around mid-rapidity is apparent in Fig. [4] at 
Tevatron energies. 
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Figure 4: Energy dependence of the average transverse momentum of 
primary charged particles in pp an d pp co llisions. Data from other 
experiments are taken from (20ll2lll2al2ll2ll23l . 

Indeed, a decrease of {pj} by about 2% is found 
between \r/\ < 0.2 and 0.6 < \rj\ < 0.8 in a pseudora- 
pidity dependent analysis of the present data. A consis- 



tent decrease of {pj) is also observed in the CMS data, 
when pseudorapidity is increased |2(| 2^]. Likewise, a 
decrease of {pj) by about 5% between \rj\ < 0.8 and 
< 2.5 is found at = 900 GeV in PYTHIA. 
Charged particle transverse momentum distribu- 
tions can be used to tune Monte Carlo event gener- 
ators of hadron-hadron interactions, such as PYTHIA 
and PHOJET. Recently, PYTHIA was tuned to describe 
the energy dependence of existing measurements, e.g. 
with respect to the treatment of multiple parton inter- 
actions and divergencies of the 2— >2 parton scattering 
cross-section at small momentum transfers. 

cv. 10f 




P T (GeV/c) 



Figure 5: Top: Comparison of the primary charged particle differen- 
tial yield in INEL pp collisions at V? = 900 GeV < 0.8) to results 
from PHOJET and PYTHIA tunes 109 (Hi, 306 Q and 320 (H. 
Bottom: Ratio between the Monte Carlo simulation and the data. The 
shaded area indicates the statistical and systematic errors of the AL- 
ICE data added in quadrature. 

In Fig. |5] the results for INEL events are com- 
pared to PHOJET and different tunes of PYTHIA, D6T 
(tune 109) 111, PerugiaO (tune 320) and ATLAS- 
CSC (tune 306) 12811 . The best agreement is found 
with the PerugiaO tune, which gives a fair description 
of the spectral shape, but is approximately 20% below 
the data. The D6T tune is similar to PerugiaO below 
2 GeV/c but underestimates the data more significantly 
at high p T . PHOJET and the PYTHIA ATLAS-CSC 
tune fail to reproduce the spectral shape of the data. 
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They overestimate the yield below 0.7 GeV/c and fall 
short of the data at high p T . We note that PHOJET 
and ATLAS-CSC agree best with the charged particle 
multiplicity distributions at t/s — 0.9, 2.36 and 7 TeV, 
respectively 

Fig.|6]shows the pj spectra in 1NEL events for three 
different multiplicity selections (n acc ) along with fits to 
the modified Hagedorn function (Eq. [TJ. A consider- 
able flattening of the tails of the spectra is visible with 
increasing multiplicity. The fit parameters pj$ and b 
drop by more than 50% from the lowest to the highest 
multiplicities. The results for the fit parameters in bins 
of n acc are listed in Table [2] along with the average mul- 
tiplicity (« c h) assigned to each « acc as determined from 
Monte Carlo simulations. 
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Figure 6: Upper panel: Transverse momentum spectra of primary 
charged particles in INEL pp collisions at yfs = 900 GeV (\rf\ < 0.8), 
normalized to the total number of INEL events /V C vt, for three dif- 
ferent event multiplicities together with the modified Hagedorn fits 
(Eq. described in the text. The fits are performed in the range 
0. 15 < pj < 4 GeV/c and extrapolated to pj = 0. The error bars in- 
dicate the statistical and systematic errors added in quadrature. Lower 
panel: Ratios of the pj spectra in different multiplicity ranges to the 
inclusive pj spectrum in INEL events. 

Also shown in Fig.|6]are ratios of pj spectra in dif- 
ferent multiplicity regions over the inclusive pj spec- 
trum in INEL events. A very pronounced multiplicity 
dependence of the spectral shape is manifest, exhibiting 
enhanced particle production at high pj in high multi- 
plicity events. At pj < 0.8 GeV/c the trend is opposite, 



albeit with a much weaker multiplicity dependence. The 
evolution of the spectral shape with multiplicity may 
shed light on different particle production mechanisms 
in pp collisions. A qualitatively similar evolution of the 
pj spectra with multiplicity has been seen in pp data at 
V^ = 200 GeV IH. 

The average transverse momentum (pj) as a func- 
tion of the multiplicity of accepted particles (rc acc ) in 
INEL pp collisions at = 900 GeV is shown in the 
left panel of Fig. [7j For all three selected pj ranges 
a significant increase of (pj) with multiplicity is ob- 
served. Most significantly for 0.5 < pj < 4 GeV/c, 
the slope changes at intermediate multiplicities. 

In the right panel of Fig. [7] the same data is shown 
as a function of « c h after application of the weighting 
procedure (Eq.0. In comparison to model calculations, 
good agreement with the data for 0.5 < pj < 4 GeV/c is 
found only for the PYTHIA PerugiaO tune (Fig. [8] left 
panel). In a wider pseudorapidity interval Qt]\ < 2.5), 
similar agreement of the data with PerugiaO was re- 
ported by ATLAS 03 • For 0.15 < p T < 4 GeV/c, Pe- 
rugiaO and PHOJET are the closest to the data, as shown 
in the right panel of Fig.[8j however, none of the models 
gives a good description of the entire measurements. 

6. Conclusion 

A measurement is presented of the primary charged 
particle transverse momentum spectrum and of the 
mean transverse momentum in pp collisions at = 
900 GeV with the ALICE detector at the LHC. Good 
agreement with previous results from LHC is found up 
to pj — 1 GeV/c. At higher pj, the data exhibit a harder 
momentum spectrum of primary charged particles than 
other measurements in pp and pp collisions at the same 
energy. We argue that this is most likely related to the 
different pseudorapidity intervals studied. The average 
transverse momentum in \q\ < 0.8 is (pj) inel = 0.483 + 
0.001 (stat.) +0.007 (syst.) GeV/c and (p T ) NSD = 
0.489 + 0.001 (stat.) ±0.007 (syst.) GeV/c. None of the 
models and tunes investigated simultaneously describes 
the pj spectrum and the correlation between {pj) and 
n c h. In particular in the low pj region, where the bulk 
of the particles are produced, the models require further 
tuning. These measurements will help to improve the 
phenomenological description of soft QCD processes 
and the interplay between soft and hard QCD. The pre- 
sented data demonstrate the excellent performance of 
the ALICE detector for momentum measurement and 
will be used as a baseline for measurements at higher 
LHC energies and for comparison with particle produc- 
tion in heavy-ion collisions. 
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Figure 7: The average transverse momentum of charged particles in INEL pp events at = 900 GeV for three different pj ranges as a function 
of 

Wacc (left panel) and as a function of (right panel). The error bars and shaded areas indicate the statistical and systematic errors, respectively. 
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Figure 8: The average transverse momentum of charged particles for 0.5 < pj < 4 GeV/c (left panel) and 0.15 < pj < 4 GeV/c (right panel) in 
INEL pp events at V* = 900 GeV as a function of n c h in comparison to models. The error bars and the shaded area indicate the statistical and 
systematic errors of the data, respectively. In the lower panels, the ratio Monte Carlo over data is shown. The shaded areas indicate the statistical 
and systematic uncertainty of the data, added in quadrature. 
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Table 1: Contributions to the systematic uncertainties on the differential primary charged particle yields l/N^ \H2npj) d 2 N c h/(di] Apr) and the 
average transverse momentum (pj) . Ranges are given if the contributions are pj dependent. 





1 1 d 2 /V ch 

yVcvi 27Tpr dl] dpr 


(Pt) 


p T range (GeV/c) 


0.15-10 


0.5-4 


0.15-4 


0-4 (extrap.) 


Track selection cuts 


0.2-4% 


negl. 


0.3% 


0.5% 


Contribution of diffraction (INEL) 


0.9-1% 


negl. 


negl. 


negl. 


Contribution of diffraction (NSD) 


2.8-3.9% 








Event generator dependence (INEL) 


2.5% 


negl. 


negl. 


negl. 


Event generator dependence (NSD) 


0.5% 








Particle composition 


1-2% 


0.1% 


negl. 


0.1% 


Secondary particle rejection 


0.2-1.5% 


negl. 


0.1% 


0.2% 


Detector misalignement 


negl. 


negl. 


negl. 


negl. 


ITS efficiency 


0-1.6% 


negl. 


0.3% 


0.5% 


TPC efficiency 


0.8-4.5% 


negl. 


0.5% 


0.7% 


SPD triggering efficiency 


negl. 


negl. 


negl. 


negl. 


VZERO triggering efficiency (INEL) 


negl. 


negl. 


negl. 


negl. 


VZERO triggering efficiency (NSD) 


0.2% 








Beam-gas events 


negl. 


negl. 


negl. 


negl. 


Pile-up events 


negl. 


negl. 


negl. 


negl. 


Total (INEL) 


3.0-7.1% 


0.1% 


0.7% 


1.0% 


Total (NSD) 


3.5-7.2% 








R weighting procedure 




3.0% 


3.0% 


3.0% 


Extrapolation to pj = 








1.0% 


Total 




3.0% 


3.1% 


3.3% 
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Table 2: Parameters of the modified Hagedorn fits (Eq.[T} to the transverse momentum spectra. The fit range in pr is 0.15 - 10 GeV/c for the 
multiplicity integrated spectra (first two rows) and 0.15-4 GeV/c for the spectra binned in multiplicity. The errors are statistical and systematic 
added in quadrature. Also given are the average multiplicites (n c h) of events contributing to the « acc bins, as determined from Monte Carlo. 



event class 


^acc 


<«ch> 


Ptx) (GeV/c) 


b 


INEL 
NSD 


all 
all 




1.05 ±0.05 
1.05 ±0.05 


7.92 ± 0.04 
7.84 ± 0.04 


INEL 


1 


2.1 ±0.1 


2 


.64 


± 





.29 


16.50 ± 1.32 


INEL 


2 


3.5 ±0.1 


1 


.86 


± 





15 


12.58 ±0.69 


INEL 


3 


4.8 ±0.1 


1 


.49 


± 





11 


10.56 ± 0.45 


INEL 


4 


6.1 ±0.1 


1 


.26 


± 





.08 


9.28 ±0.34 


INEL 


5 


7.4 ±0.1 


1 


.16 


± 





.07 


8.60 ±0.28 


INEL 


6 


8.7 ±0.1 


1 


.04 


± 





.06 
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